Sympathetic cooling of polyatomic molecules with S'-state atoms in a magnetic trap 
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We present a rigorous theoretical study of low-temperature collisions of polyatomic molecular 
radicals with So atoms in the presence of an external magnetic field. Accurate quantum scattering 
calculations based on ab initio and scaled interaction potentials show that collision-induced spin 
relaxation of the prototypical organic molecule CH2(X 3 -Bi) (methylene) and 9 other triatomic rad- 
icals in cold 3 He gas occurs at a slow rate, demonstrating that cryogenic buffer-gas cooling and 
magnetic trapping of these molecules is feasible with current technology. Our calculations further 
suggest that it may be possible to create ultracold gases of polyatomic molecules by sympathetic 
cooling with alkaline-earth atoms in a magnetic trap. 
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The study of molecular structure and interactions at 
low temperatures is a promising route toward exploring 
novel phenomena in physics and chemistry, with applica- 
tions ranging from quantum simulation of strongly cor- 
related many-body systems [l| to quantum information 
processing [i[ and metrology Q. Recent experimental 
observation of chemical reactions in ultracold molecular 
gases Q and atom-molecule mixtures [B| has opened a 
new chapter in chemical physics by demonstrating the 
possibility of controlling bimolecular chemical reactions 
with external electromagnetic fields Q. Some of the 
most important experimental techniques for the creation 
of cold molecules include photoassociation of ultracold 
atoms [H, IH , Stark deceleration of molecular beams Q , 
and cryogenic buffer gas cooling |8. 

With few exceptions [H,[lo|], previous experimental and 
theoretical work has been limited to diatomic molecules 
[H . Cooling and trapping of polyatomic molecules would 
greatly enrich the scope of molecular physics and allow 
the study of complex chemical reactions |12||. molecular 
decoherence [l3[ and parity violation |l4| in the well- 
controlled environment of an electromagnetic trap. A 
potentially viable technique for creating ultracold poly- 
atomic molecules could be based on cryogenic buffer gas 
cooling Q or Stark deceleration [l(| followed by magnetic 
trapping and subsequent sympathetic cooling with ultra- 
cold atoms Low-temperature collisional prop- 
erties of polyatomic molecules play a key role in this 
scheme: only the molecules in low-field-seeking Zeeman 
states can be confined in a magnetic trap, and such 
molecules are vulnerable to spin relaxation induced by 
collisions with buffer gas atoms, causing trap loss. For 
efficient sympathetic cooling, the number of elastic colli- 
sions per inelastic collision 7 must exceed 100 [H, [H]. 

In this Letter, we use rigorous quantum scattering cal- 
culations to demonstrate the attainability of sympathetic 
cooling of 10 polyatomic molecular radicals with cold 3 He 
atoms in a magnetic trap. To our knowledge, this is 
the first rigorous theoretical analysis of low-temperature 



collisions of polyatomic molecular radicals that incorpo- 
rates in a rigorous manner both intra- and intermolecu- 
lar interactions and the effects of an external magnetic 
field on collision dynamics. Our calculations provide 
compelling evidence that methylene (CH2) — the proto- 
typical organic molecule — can be cooled and magneti- 
cally trapped by collisions with cryogenic He buffer gas, 
thereby opening up novel opportunities for research in 
astrophysics, combustion, organic chemistry, molecular 
interferometry, and precision spectroscopy. In addition, 
our study provides insight into the mechanisms of spin 
relaxation of polyatomic molecules in collisions with S- 
state atoms, demonstrating that polyatomic molecules 
with one unpaired electron undergo inelastic spin relax- 
ation much slower than molecules with two unpaired elec- 
trons. Our calculated low-temperature collision rates for 
polyatomic molecules are similar to those measured pre- 
viously for the diatomic molecules CaH( 2 £) and NH( 3 I]), 
which indicates that advanced collisional cooling tech- 
niques that are currently under development for diatomic 
molecules [l5-17| can be extended in a straightforward 
manner to polyatomic molecules, thereby demonstrating 
a pathway to ultracold ensembles of chemically diverse 
molecules that cannot be created by any other cooling 
method in existence today. 

We begin by specifying the Hamiltonian for an open- 
shell polyatomic molecule in the ground electronic and 
vibrational states colliding with a 1 Sq atom in the pres- 
ence of an external magnetic field (H = 1) 
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where \i is the reduced mass of the collision complex, R 
is the atom-molecule separation vector of length R, £ 2 is 
the orbital angular momentum for the collision and -ff mo i 
describes the internal structure of the molecule 

H mol = AN 2 + CN 2 + BN 2 + H cd + H SI + H ss + H cxt (2) 

where A. B. and C are the rotational constants of an 
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asymmetric top, N a (a = x, y, z) yield the components 
of rotational angular momentum N in the molecule-fixed 
(MF) coordinate frame. We identify the MF axes with 
the principal axes of inertia of the molecule [l^- The 
spin-rotation interaction is given by [l9| 



H sr =^N-S+J2 

q=-2 



h x -iy)[v 2 q2 {h) + vl_ 2 {ii)} 



[N®S]& (3) 



where S is the electron spin, j a are the spin-rotation 
constants, 7 = \{"f x + l y + 7a) j T^qq'i^) are the Wigner 
D-functions, and ft are the Euler angles which specify the 
orientation of MF axes in the space-fixed (SF) coordinate 
frame. The spin- spin interaction in S = 1 molecules may 
be written as [19l ulfl 



H ss = Y 

q=-2 



-{E + D)[V 2 q2 (h)+Vl_ 2 (h)\ 
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where D and E are the zero-field-splitting parameters 
(TsL [20^ . The interaction of the molecule with an external 
magnetic field is described by H cxt = 2fJ-oBSz, where B 
is the magnetic field strength, /xo is the Bohr magneton 
and Sz yields the SF projection of S. The term H c a 
accounts for the effects of centrifugal distortion (l9| . 

To solve the quantum scattering problem for the 
Hamiltonian |T]), we expand the wave function of 
the collision complex in the fully uncoupled SF basis 
\NM N K N )\SM s )\im e ), where \NM N K N ) are the sym- 
metric top eigenfunctions, M/v, mi, and M$ are the SF 
projections of N, £, and S, and is the MF projec- 
tion of N . This expansion leads to a system of coupled- 
channel (CC) equations parametrized by the matrix el- 
ements of the Hamiltonian (fT]), which can be evaluated 
analytically [2l| using the known spectroscopic constants 
of polyatomic radicals |igL [2l| . The matrix elements of 
the atom-molecule interaction potential can be evaluated 
by expanding the potential in angular functions [18l 121 1 
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YPi x ^) + (-rV^^(Cl)]Y Xmx (R) (5) 



where R = R/R. 

At this point, we specialize our analysis to the 
He + CH2 collision system. CH2 (methylene) is the sim- 
plest polyatomic molecular radical with a triplet ground 
state of X 3 Bi symmetry, which serves as a prototype 




Fig. 1. A contour plot of the ab initio He-CH2 PES in Carte- 
sian coordinates Xmf = R sin 6 cos cj> and Zmf = R cos 6 for 
4> = 90°. Energies are in cm" 1 (1 cm" 1 = 1.439 K). 



for studying the mechanisms of inelastic collisions and 
chemical reactions of polyatomic molecules [24[. Methy- 
lene has recently been detected in interstellar space [25|, 
and in Titan's atmosphere, where it is formed by pho- 
todissociation of methane [26(. In addition, CH2 is the 
simplest known carbene, an important class of reaction 
intermediates in organic chemistry [27j. 

The ground rotational state of CH2 is split by an ex- 
ternal magnetic field into three Zeeman sublevels with 
A/5 = and ±1. Like homonuclear molecules with non- 
zero nuclear spins, methylene can exist in two spin mod- 
ifications (ortho and para), only one of which (0-CH2) 
is amenable to magnetic trapping (2l| . In this work, we 
consider collisions of 0-CH2 molecules in their maximally 
spin-stretched Zeeman sublevel of the ground rotational 
state I Ms = 1). Collisions with He atoms induce spin re- 
laxation transitions to the Zeeman states \Ms = 0) and 
I Ms = — 1), which are the main focus of this work. To 
quantify the dynamics of spin relaxation, we solve the 
CC equations numerically on a grid of R from 2 to 40ao 
with a grid step of 0.04 ao, The CC basis set included 
five rotational states of CH2, and six partial wave states 
(£ = - 5), leading to 1123 CC equations for M = 0. 
The calculated cross sections were converged to <10%. 

To evaluate the potential energy surface (PES) for 
the He-CH2 collision complex, we use a high-level ab 
initio approach based on the coupled-cluster method 
with single, double, and non-iterative triple excitations 
[CCSD(T)] [lU and a large basis set of cc-pVTZ qual- 
ity [23|]. A total of 985 ab initio points were calculated 
and fitted by an analytic function consisting of six sets of 
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Fig. 2. (a) Cross sections for elastic scattering (dashed line) 
and spin relaxation in He + CH2 collisions plotted vs collision 
energy for different magnetic field strengths: O.Of T (circles), 
0.1 T (diamonds), and 1 T (triangles). The s-wave scattering 
lengths for B = 0.1 T are 18.6 (f s = 1) and 13.9 a {f s = 10). 



long-range pair potentials. The root-mean-square error 
of the fit was 0.8 cm -1 . A contour plot of the calcu- 
lated PES in shown in Fig. [TJ The global minimum is 
located at R = 7.14 ao and 9 = it and has a well depth 
of 6.3 cm . The energy difference between the linear 
and T-shaped configurations amounts to ^1.5 cm -1 for 
4> = 90°, demonstrating that the He-CPi2 interaction is 
weakly anisotropic. 

Figure fj] shows the cross sections for elastic scatter- 
ing and spin relaxation in He + CH2 collisions as func- 
tions of collision energy. A pronounced maximum near 
Ec ~ 0.1 K occurs due to a shape resonance in the inci- 
dent collision channel. At ultralow collision energies, the 
cross sections for elastic energy transfer approach a con- 

— 1/2 

stant value and those for spin relaxation vary as E c . 
according to the Wigner threshold law for s-wave scat- 
tering. The sensitivity of the cross sections to magnetic 
field is the strongest at ultralow collision energies, and 
near the shape resonance at Ec ~ 0.1 K. 

Figure H^a) shows the calculated ratio of the rate 
constants for elastic scattering and spin relaxation in 
He + CH2 collisions at a magnetic field of 0.1 T. The 
ratio remains high (7 > 10 4 ) over a wide range of tem- 
peratures from 1 K down to 10 -5 K. In order to ex- 
amine the sensitivity of these results to the anisotropic 
part of the interaction potential, we repeated scattering 
calculations with a modified He-CH2 interaction poten- 
tial obtained by scaling the V\^(R) terms with A > 
(Eq. |5} by a constant factor f s = 10. The results 
plotted in Fig. \3j[&) demonstrate that even for this sig- 
nificantly more anisotropic interaction potential, j(T) 
hardly ever falls below 100, demonstrating that CH2 
molecules trapped in the presence of cold 3 He gas will 
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Fig. 3. Ratios of elastic to inelastic collision rates 7(T) for (a) 
He + CH2 and (b) He + NH2 as functions of temperature. 
Full lines: results for the unmodified He-CH2 interaction po- 
tential (f s — 1), dashed lines: results for the scaled potential 
[H. The magnetic field is 0.1 T. 



be extremely stable against collisional spin relaxation. 

In order to explore the possibility of sympathetic cool- 
ing and magnetic trapping of polyatomic molecules other 
than CH 2 , we have extended our analysis to include the 
NH 2 molecule in its ground rovibronic state of 2 A sym- 
metry. Figure E3h) shows the elastic-to- inelastic ratio for 
He + NH2 collisions computed using the He-CH2 interac- 
tion potential and the accurate spectroscopic constants of 
NH 2 " 



2l| . The interactions of S-state molecules with He 



atoms are generally weakly anisotropic [H, [Hi , so we be- 
lieve that the replacement of the unknown He-NH2 inter- 
action PES by the He-CH2 PES computed in this work is 
a reasonable approximation. The elastic-to-inelastic ra- 
tio for He + NH2 is remarkably high, exceeding that for 
He + CH2 by nearly 7 orders of magnitude. Calculations 
with the modified PES yield 7 ~ 10 4 , still perfectly ac- 
ceptable in the context of buffer gas cooling and magnetic 
trapping experiments 

Table I lists the values of 7 for 18 molecular radicals 
obtained from quantum scatte ring calculations using the 
accurate molecular constants [2l| and both the original 
and modified He-CH2 PESs. Because the magnitude of 
7 can be sensitive to the anisotropic part of the PES 
111, the f s = l results are best thought of as order-of- 
magnitude estimates, whereas the f s > 10 results may 
be regarded as conservative lower bounds [2j| . A few im- 
portant qualitative conclusions can be drawn based on 
the results listed in Table I. All S — 1/2 molecules ex- 
cept HCO have extremely large elastic-to-inelastic ratios, 
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TABLE I. Ratios of the rate constants for elastic energy trans- 
fer and spin relaxation calculated for 18 polyatomic molecu- 



lar radicals with He using unmodified (f a 



and strongly 



anisotropic (f s — 10 unless noted otherwise [29() interaction 
potentials for He-CH 2 . T = 0.5 K and B = 0.1 T. 



Molecule 


7 




Molecule 


7 




(S - 11 


f — 1 


f — in 

J s — 


(S - i"l 


f — 1 


f — 10 

Js — 1U 


CH 2 


4.5 x 10 4 


326 


NH 2 


1.1 x 10 12 


3.6 x 10 4 


CHF 


1.7 x 10 3 


186 


PH 2 


4.2 x 10 10 


6.3 x 10 3 


CHC1 


418 


11 


AsH 2 


9.5 x 10 9 


140 


CHBr 


156 


3 


H0 2 


1.8 x 10 9 


5.6 x 10 3 


CHI 


64 


10 


HCO 


1.9 x 10 11 


25 


CF 2 


80 


11 


NF 2 


4.6 x 10 7 


894 


CC1 2 


11 


9 


N0 2 


2.1 x 10 8 


541 


CFC1 


10 


8 


C10 2 


4.9 x 10 4 


260 


SiH 2 


2.5 x 10 4 


34 


CH 3 


9.8 x 10 12 


2.8 x 10 7 



which makes them amenable to sympathetic cooling and 
magnetic trapping experiments using cold 3 He gas. This 
conclusion is independent of whether unmodified (/ s = 1) 
or strongly anisotropic (f s > 10) interaction potentials 
are used in scattering calculations, even though the mag- 
nitude of 7 is sensitive to f s . Apart from CH2 , only one 
triplet radical (CHF) satisfies the criterion 7 > 100 for 
f a = 10. To explain this result, we note that spin re- 
laxation in S = 1 molecules is driven by the spin-spin 
interaction (0)), which couples the ground initial Zeeman 
state \N = 0, A/j = 1) to a manifold of rotationally ex- 
cited states \N = 2, Ms). These couplings lead to di- 
rect spin-flipping transitions, whose cross sections scale 
quadratically with D [ll[. In the absence of the spin- 
spin interaction, spin relaxation occurs via a much less 
efficient indirect mechanism [28| . 

In summary, we have presented a rigorous theoreti- 
cal analysis of low-temperature collisions of polyatomic 
molecular radicals with closed-shell atoms in the presence 
of an external magnetic field using an accurate ab initio 
interaction potential. The calculations demonstrate that 
spin relaxation in 3 He + CH2 collisions occurs at a re- 
markably slow rate of 1.2 x 10~ 14 cm 3 /s at T = 0.5 K, 
thereby making CH2(X 3 £>i) an ideal candidate for sym- 
pathetic cooling experiments using cold 3 He gas. We 
have also presented model calculations of spin relaxation 
rates for other polyatomic molecules, demonstrating that 
two S = 1 molecules and eight S =1/2 molecules listed 
in Table I should be amenable to cryogenic buffer-gas 
cooling and magnetic trapping with long lifetimes. 

This work demonstrates the feasibility of sympathetic 
cooling and magnetic trapping of a large class of molec- 
ular radicals of importance in organic chemistry, astro- 
physics, precision spectroscopy, and molecular interfer- 



ometry, which may open up novel avenues of research. 
Inelastic collisions, chemical reactions and three-body re- 
combination of polyatomic molecules can now be stud- 
ied in cold collision experiments @, Accurate mea- 
surements of radiative lifetimes and reaction rates of co- 
trapped molecular radicals would greatly facilitate as- 
trochemical modeling of dense cold interstellar clouds 
[l2l . [25| . Magnetic trapping of polyatomic molecules 
may pro vide novel routes to probing quantum decoher- 
ence 



13( and controlling reaction mechanisms with ex- 
ternal fields [fj. Finally, we note that the calculated 
elastic-to-inelastic ratios for molecular species listed in 
Table I are similar in magnitude to those already mea- 
sured for CaH( 2 E) [H] and NH( 3 E) OH, which indicates 
that it may be possible to create ultracold ensembles of 
polyatomic molecules via sympathetic cooling with laser- 
cooled alkaline earth or spin-polarized N atoms 
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